arXiv:1508.01005vl [cond-mat.mtrl-sci] 5 Aug 2015 


Atomistic details of oxide surfaces and surface oxidation: the example 

of copper and its oxides 

Chiara Gattinoni and Angelos Michaelides 

Thomas Young Centre, London Centre for Nanotechnology and Department of Chemistry 

University College London 
17-19 Gordon Street, London, WC1H OAH, UK 


Abstract 

The oxidation and corrosion of metals are fundamental problems in materials science and tech¬ 
nology that have been studied using a large variety of experimental and computational techniques. 
Here we review some of the recent studies that have led to significant advances in our atomic-level 
understanding of copper oxide, one of the most studied and better understood metal oxides. We 
show that a good atomistic understanding of the physical characteristics of cuprous (CU 2 O) and 
cupric (CuO) oxide and of some key processes of their formation has been obtained. Indeed, the 
growth of the oxide has been proven to be epitaxial with the surface and to proceed, in most cases, 
through the formation of oxide nano-islands which, with continuous oxygen exposure, grow and 
eventually coalesce. We also show how electronic structure calculations have become increasingly 
useful in helping to characterise the structures and energetics of various Cu oxide surfaces. However 
a number of challenges remain. For example, it is not clear under which conditions the oxidation of 
copper in air at room temperature (known as native oxidation) leads to the formation of a cuprous 
oxide film only, or also of a cupric overlayer. Moreover, the atomistic details of the nucleation of the 
oxide islands are still unknown. We close our review with a perspective on future work and discuss 
how recent advances in experimental techniques, bringing greater temporal and spatial resolution, 
along with improvements in the accuracy, realism and timescales achievable with computational 
approaches make it possible for these questions to be answered in the near future. 


1 Introduction 

Copper is a material which has accompanied human pre-history and history, and it is still highly 
relevant today. Cold working of copper has been performed for at least 10,000 years and smelting of 
copper ore for around 7,000 years 1 . Its use as construction materials, e.g. in piping, can be dated 
back to ancient Egypt, and its importance in this field has not diminished nowadays. In the modern 
world it has acquired further uses, for example in electrical systems and electronic devices. 

Within this long history, the properties of copper have been extensively studied and exploited, 
however much is still unknown about this important metal. In particular, the oxidation and corrosion 
of copper, which impacts its performance in industrial and technological applications, is still not 
completely understood. 

Copper is found to readily oxidise at room temperature 2 ^, and the presence of an oxide layer, how¬ 
ever thin, can compromise its uses in technology. As an example, copper could be an environmentally 
friendly and low-cost substitute for the (currently used) tin-lead or (promising) gold- and silver-based 
solder alloys in electronic packaging, if there was a way to prevent its oxidation 5 . Moreover, copper 
canisters are used for nuclear waste disposal, and understanding the oxidation and corrosion of copper 
in anaerobic conditions is thus really important 6 . On the other hand, the existence of stable copper 
oxides at room temperature, with a ~ 2.0 eV band gap, makes them interesting for catalytic^, gas 
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sensing®, optoelectronic and solar technologies 9 10 . Thus, there currently is a two-fold interest in un¬ 
derstanding copper oxides: from the one hand, to mitigate against technological failure, on the other 
hand, to exploit their potential industrial applications. 

Copper is considered a model system to understand the formation of metal oxides in general. 
The atomistic details of the oxidation process tend to be system-specific, with some metals showing 
uniform oxide growth ( e.g. Ref. 1111 ), other complex temperature-dependent phenomena such as surface 
roughening (e.g. Ref. 12 ) and island formation (e.g. Ref. 13235 ). However, the copper oxidation process 
is one of the most studied with a large number of experimental and computational methods, and one 
of the better understood. Therefore, a detailed understanding of copper oxidation, of the techniques 
used to study it and of the challenges which are still open is invaluable when considering the oxidation 
process on any other system. 

In this review, we discuss the status of knowledge of copper oxidation from the atomistic point 
of view, which we believe is of key importance if we want to learn how to prevent or manipulate 
copper oxidation. We cannot hope to provide a complete review of all the work done on this subject 
since the beginning of the last century 1 ® 16 . We are therefore only going to focus on recent surface 
science, spectroscopy and atomistic computational work which has been performed to understand the 
properties of copper oxides and their formation, and on the open challenges that can be addressed 
using these techniques. 

First, a brief overview on the experimental and computational techniques which have been used 
for the study of oxide structures and oxidation kinetics is given (Sec. [ 2 ]), in order to clarify some of 
the terminology used throughout the review. The structural and electronic characteristics of the bulk 
oxides and their surfaces are then presented (Sec. [3]). We subsequently look at the interaction of clean 
copper surfaces with oxygen and the initial stages of controlled oxidation (Sec. |4j) as well as long-term 
oxidation (Sec. [5]). Finally, in Sec.[6j conclusions and perspectives are given. 

We hope that it will be clear from the following that tremendous progress has been made in 
understanding the atomistic details of copper oxides and their formation under different conditions. 
However, equally important gaps in our understanding remain, especially in terms of the formation 
kinetics and the structure of the resulting oxide surfaces. 

2 Experimental and computational techniques 

An enormous number of experimental 1 ® 20 and computational techniques are available to investigate 
the physical and chemical characteristics of solids, surfaces and surface kinetic processes. Many of 
them have been used over the years to e.g., understand oxide structures, characterize oxide surfaces, 
understand the oxidation kinetics and investigate bulk properties of the oxides. For clarity, in this 
section we provide a brief introduction to the most relevant techniques which have been used on 
copper oxide and which are going to be mentioned in the following sections, with an emphasis on 
strengths and weaknesses of each approach. 

In early studies of oxidation, thermogravimetric analysis (TGA), where changes in physical and 
chemical properties of materials are measured as a function of time, has been widely used to study 
the onset of oxidation by recording the mass gain of a sample under oxygen exposure. Whilst useful 
in providing a broad overview of the extent of oxidation, this technique is however unable to provide 
atom-resolved information. 

Imaging techniques such as electron microscopy (EM), transmission electron microscopy (TEM) 
and their derivatives (e.g., high-resolution transmission electron microscopy (HRTEM) or field emission 
scanning electron microscopy (FESEM)) and surface-specific techniques such as scanning tunnelling 
microscopy (STM) can be applied to surfaces and provide atomistic level structural information. As 
we will see, they have been amply used to image adsorption of oxygen on the copper surface, surface 
reconstructions and initial oxide formation. The atomic composition and oxidation states of the atoms 
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in a material can be obtained using spectroscopic techniques such as X-ray diffraction (XRD), Auger 
electron spectroscopy (AES), electron energy loss spectroscopy (EELS), low-energy electron diffraction 
(LEED), reflection high-energy electron diffraction (RHEED), X-ray photoelectron spectroscopy (XPS 
or ESCA) and X-ray absorption spectroscopy (XAS). The space- and time-resolution of these techniques 
has greatly improved since they were first used in this field (in the 1980s), as well as their range of 
applicability: for example, experiments at relatively high pressures can be performed nowadays 7 . 

Growth of the oxide and its atomic composition has been extensively studied by means of ellip- 
sometry. This technique measures changes in polarization as light interacts with an object and the 
resulting data are fitted with a ‘guess’ model for the material (e.g. a two-layer CU 2 O/CU model or a 
three-layer Cu0/Cu20/Cu one). 

Computer simulations have also been widely applied to study the structural, optical and vibrational 
characteristics of copper oxides. When carrying out simulations of materials at the atomic scale, 
classical empirical potentials (force fields) or more sophisticated quantum (ab initio) approaches, such 
as density functional theory (DFT), can be used. Force fields are parametrised empirical potentials 
tuned to reproduce the interactions of the atoms in the system at hand. For certain problems force 
fields can provide a faithful description. However, the reliability and transferability of such calculations 
depends primarily on the extent and quality of the data used in their construction. Moreover, force 
fields cannot generally be applied to study chemical reactions. Current development in force fields 
are addressing these issues. Indeed, parametrisations obtained by fitting large data sets using e.g., 
neural network 22 and machine learning methods 2 ®, to name just two, are improving accuracy and 
transferability. Moreover, bond order potentials (such as ReaxFF 24 ) and polarisable force fields^- 23 'are 
making it possible to simulate chemical reactions. It is difficult to make general statements about the 
sizes of systems that can be explored with various methods. However, with many standard force fields 
it is now possible to examine systems with ~ 100,000 atoms on a routine basis. In addition, it is 
possible to explore the evolution of a system of this size, again on a routine basis, with an approach 
such as molecular dynamics for several hundreds of nanoseconds. 

More accurate approaches are ab initio methods which aim to study the structure and properties 
of a material by seeking (approximate) solutions to quantum mechanical equations such as the rnany- 
electron, many-atom Schroedinger equation. These methods are more general and do not require 
system-specific parametrisations, however they are computationally more expensive. One of the most 
widely used methods to study the properties of bulk materials and surfaces is density functional theory 
(DFT) 27 30 . With DFT it is possible to make genuine predictions about structural properties of a 
material within a few percent of the experimental value. DFT also provides access to the electronic 
structure of the systems being considered and related spectroscopic properties. In DFT, the energy of 
the electronic system is determined from the electronic density by solving Schroedinger-like equations. 
Whilst exact in principle, in practice approximations have to be introduced since the functional form of 
the electron-electron interaction, called the exchange-correlation (XC) functional, is unknown. Many 
approximate XC functionals have been developed-® 4 , the most common ones being the local density 
approximation (LDA) and the generalised gradient approximation (GGA). A number of deficiencies 
in DFT arise from these approximations and the choice of an appropriate XC functional is critical in 
order to obtain meaningful results. As an example, in strongly correlated systems like CuO (typically, 
where d and / orbitals are localised), the GGA and LDA functionals provide a poor description of 
the electronic and bulk crystal structure. Moreover, the band gap obtained using GGA or LDA in 
semiconductors or insulators is generally underestimated. In these cases it is possible to add simple 
but somewhat ad hoc corrections to the functionals (e.g. the Hubbard U®, self-interaction correction 
SIC 33 ) or to use more sophisticated hybrid functionals (e.g. HSE06^, PBEO 35 ) which incorporate 
some exact Hartree-Fock exchange. Compared to standard empirical potential methods, DFT is much 
more computationally demanding, and, on a routine basis, systems with only a few 100 atoms can be 
examined and the dynamics of such systems explored for only a few tens of picoseconds. 
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Figure 1: Ball-and-stick model of a unit cell of a) cuprous oxide CU 2 O and b) cupric oxide CuO. Red 
balls represent oxygen and brown balls copper. The unit cell is shown in green. 

Other computational methods which have been used in copper oxide simulations are the GW 
method 3 ^, used when optoelectronic properties are of interest, since it is more accurate at predicting 
band structures than standard DFT XC functionals. In addition there have been a number of Hartree- 
Fock® studies and calculations with simpler approaches, such as the linear combination of atomic 
orbitals (LCAO) and tight binding methods 38 . 


3 Oxide structures 

We are going to introduce here the bulk and surface structures of the two main copper oxides, CU 2 O and 
CuO. Knowledge of the physical properties of these materials, and especially of the surface structures, 
is relevant background when trying to understand the formation and growth of the copper oxide. 

The two most common forms of the oxide (shown in Fig. [T]) are cuprite (or cuprous oxide, CU 2 O), the 
principal oxide at low temperatures and pressure, and tenorite (or cupric oxide, CuO), dominant at high 
temperatures and pressures® . Another copper oxide structure, paramelaconite (CU 4 O 3 ), exists as a 
rare mineral found in hydrothermal deposits of copper. Cuprite has long been known to be the primary 
oxide for copper at ambient conditions and there is considerable interest in its application to catalysis, 
optoelectronics and gas sensing and therefore a large amount of work has been done to determine 
its physical and chemical characteristics. Tenorite has been studied less and still relatively little is 
understood about the structure of its surfaces, with only a handful of experimental and computational 
studies performed to that aim. 

3.1 Cu 2 0 bulk properties 

In cuprous oxide (CU 2 O, see Fig. [lja,) , a cubic crystalline solid, copper has a Cu 1+ oxidation state. It is 
a p-type semiconductor with a direct band gap of 2.02 — 2.17 eV and an optical gap of 2.62 eV 10 40 41 . 
It is a promising material for a variety of industrial applications because of its band gap and because 
it shows negative thermal expansion 4243 . 

The properties of cuprous oxide have been extensively studied using empirical potentials 4 ^, tight 
binding 45 and ab initio methods. Ab initio studies include the use of DFT methods (see Table[l]), and 
good agreement with the experimental bulk structure (i.e. no more than 2.5% discrepancy between the 
calculated and experimental value of the lattice constant) and vibrational modes 46 ^ have been found. 
However, the band gap is underestimated with standard DFT XC functionals, yielding values between 
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Figure 2: Structure of the Cu-termiuated Cu20(100) and O-terminated Cu20(lll) surfaces. 0 atoms 
are shown in red and Cu atoms in brown. Atoms below the surface layer are depicted in grey in the top 
views and enclosed in a box in the side views, a) Top and side view of a copper terminated Cu 2 0(100) 
slab, consisting of alternating layers of O and Cu atoms. The surface copper atoms ‘sink’ towards 
the O-atoms plane and form Cu-Cu dimers (indicated by dl). The second dimer formed in order to 
obtain the (3\/2 x \/2)i245° reconstruction is labelled d2. b) Top and side view of the O-terminated 
Cu 2 0(lll) surface. One O-Cu-O trilayer is shown between the dashed lines in the side view. The 
surface layer presents a hexagonal structure with an unsaturated copper atom Cucus i n the middle of 
each hexagon. 


0.5 and 0.8 eV 47-49 and grossly overestimated with Hartree-Fock methods at 9.7 eV 50 . Corrections to 
DFT ( e.g . DFT+U or the GW approximation applied to DFT-GGA) give rise to band gaps that are 
in better agreement with experiment 51-5 ®, as summarised in Table [l} 

3.2 Cu 2 0 surfaces 

Experimental data are available only for Cu2O(100) and Cu20(lll) low index surfaces of cuprous 
oxide, by means of XPS (which monitors the Cu to O ratio on the surface), LEED in ultra high 
vacuum (which identifies the periodicity of the surface), and in the case of Cu20(lll), also with TEM. 

We first consider the polar Cu2O(100) surface, which is formed by a succession of purely Cu or O 
layers and yields bulk-truncated surfaces which are either O- or Cu-terminated (Fig. [2^,). Polar surfaces 
generally reconstruct 69 and three different reconstructions were observed on Cu 2 0(100) according to 
the preparation method® 0 . The most stable surface in vacuum is a Cu-terminated surface with a 
(3^2 x ^/2)RA5° reconstruction. The LEED pattern (which however showed many missing spots 
and the reconstruction is presented as tentative) was interpreted as a relaxation of the top layer of 
copper cations forming two surface dimers (dl and d2 in Fig. [2^,). Two transient structures were 
also identified: at 900 K, a (y/2 x \Z2)R45° reconstruction presenting a 1/2 layer of terminal oxygen, 
and, after long oxygen exposures, the stoichiometric (lxl) O-terminated surface. Both structures 
revert to the (3\/2 x y/2)R45° reconstruction upon further annealing over 500 K. The (3\/2 x V2)RA5° 
reconstruction has been reproduced computationally®®, and it is more stable than the Cu- or O- 
terminated (y/2 x v / 2)7?45° reconstruction 63 72 73 . 

The other experimentally studied surface is Cu20(lll). A Cu20(lll) slab is formed by a succession 
of O-Cu-O trilayers, ending with either a (stoichiometric, non polar) oxygen or (non-stoichiometric) 
copper termination (Fig. [2jo) . Photoemission and LEED experiments on this Cu-terminated (111) 
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Ref. 

a [A] 

B [GPa] 

Band gap 
[eV] 

XC 

Cortona^ 

4.221 

141 

- 

LDA 

Filippetti^ 1 

4.23 

- 

0.55 

LDA 

Gordienko^ 1 

4.2696 

- 

2.87 

LDA 

NitP 

4.216 

- 

0.52 

LDA 

Heinemann^ 

4.1656 

- 

0.99 

LDA+U (U-J=6.52) 

Tran^l 

4.27* 

- 

0.63-0.94 

LDA+U (3 <U-J< 11) 

Cortona 1 ^ 

4.359 

106 

- 

PBE 

Islam^ 

4.312 

- 

0.7 

PBE 

Isseroff^ 

4.18 

145 

0.68 

LDA 

Isseroff^ 

4.10-4.17 

135-143 

0.81-1.15 

LDA+U (2 <U-J< 8) 

Isseroff^ 

4.31 

109 

0.43 

GGA 

Isserofl 1 ^ 1 

4.26-4.30 

96-106 

0.54-0.84 

GGA+U (2 <U-J< 8) 

Isseroff^ 1 

4.28 

114 

2.84 

PBE0 

Isserof?^ 1 

4.29 

114 

2.04 

HSE 


4.317 

- 

- 

PBE 

Martinez- 

Ruizf 48 ! 

4.3 

108 

0.5 

PBE 

Bohnen 1 ^ 1 

4.30 

112 

- 

PBE 

Soon— 

4.34 

104 

- 

PBE 

Soon 1 ^ 1 

4.32 

104 

0.64 

PBE 

RuizP^ 

4.435 

100 

9.7 

HF 

RuizfMl 

4.277 

93 

- 

HF+LYP 

Heinemann^ 

4.2675 

- 

2.02 

HSE06 

Scanlon*^ 

- 

- 

2.12 

HSE 

Tran^l 

4.27* 

- 

0.79-2.77 

PBE0 

Bruneval*^ 

- 

- 

1.34 

G 0 Wo 

Bruneval^ 

- 

- 

1.97 

scGW 

Filippetti^ 1 

4.23 

- 

1.8 

SIC 

Kotani^ 

- 

- 

1.97 

scGW 

Lany 1 ^ 

- 

- 

2.03 

GW+V d 

Experiment ll MF 

4.27 

106-138 

2.02-2.17 

- 


Table 1: Summary of the calculated lattice constant a, bulk modulus B and band gap for CU 2 O. Infor¬ 
mation on the exchange-correlation (XC) functional is also given. PBE is a type of GGA functional. 
scGW, GOWO and GW+V^ are different applications of the GW method. HF+LYP is a Hartree-Fock 
type calculation with an a posteriori correction using the LYP correlation functional. For calculations 
using the +U Hubbard correction, values of the relevant parameters used (U-J) are given ion eV units. 
*These calculations were performed at the experimental value of the lattice constant. 
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surface found, after annealing in UHV, a nearly stoichiometric reconstruction, with (\/3 x \/3)i?30 o 
periodicity, attributed to an ordered 1/3 of an atomic layer of oxygen vacancied 54 76 -. Conversely, 
annealing in oxygen gives a stoichiometric oxygen-terminated surface (with possibly the loss of the 
unsaturated copper atoms, Cu cus hi Fig. [ 2 ) 3 , at the centre of the hexagons/ 54 . The stability of both 
polar and non-polar stoichiometric Cu 20 (lll) surfaces has been studied also computationally 47 74 78 
by means of DFT combined with, in some cases, ab initio thermodynamics. Ab initio thermodynam¬ 
ics (see e.g. 79 80 ) is a technique which allows one to estimate relative system stabilities in different 
environments; in this case, the stability of different oxide surfaces at a range of temperatures and 
oxygen partial pressures. These studies showed that the experimentally observed Cu20(lll) surfaces 
are indeed the most stableJ 76 . 

A wider range of low-index surfaces has been probed using DFT methods than experimentally' 7 . 
The predicted lowest energy structure in high O 2 pressure is the already mentioned Cu20(lll) surface 
missing the unsaturated Cu cus atom and at low O 2 pressure the Cu2O(110) surface with a CuO-like 
surface reconstruction, which has not been observed experimentally (the (110) surface has not been 
studied yet). 

It is thus evident that there is scope for further experimental work on this topic, to confirm the 
computational predictions or to propose new reconstructions. Surface-sensitive techniques such as STM 
and LEED, which allow a direct imaging of the surfaces, would be best suited to this aim. Moreover, 
theoretical methods would be well suited to understand the reason behind the transitions between 
different reconstructions at different conditions of temperature and pressure. 

3.3 CuO bulk properties 

In CuO, the copper atom has oxidation state Cu 2+ . The unit cell has monoclinic symmetry® (see 
Fig. 0 and it contains four CuO dimers in the crystallographic unit cell, and two CuO units in the 
primitive cell. Each copper atom is located in the centre of an oxygen parallelogram. Each oxygen 
atom, in turn, has a distorted tetrahedral copper coordination. CuO is a p-type semiconductor- 2 ’ 185 
and it is antiferromagnetic in its ground state 86 89 . 

From a DFT point of view, standard XC functionals alone are not accurate enough to reproduce the 
distorted nature of the lattice, and, upon structural relaxation, they produce an orthorhombic rather 
than a monoclinic structure 90 . However, the use of more sophisticated functionals 52 60 90 93 allows for 
the reproduction of the triclinic structure and good agreement with experiments on structural and 
vibrational data®H 9 ^ (see Table [ 2 ]) . 

The research focus for CuO has been for a long time on its optical properties, for phototermal 
and photoconductive applications, and on magnetic and phase stability properties, for possible high- 
temperature superconducting applications. The electronic properties of CuO, instead, have been stud¬ 
ied only more recently. The band gap of CuO has been experimentally measured to be 1.2 — 1.9 
e Yl 10 l 41 l 84 l 85 l 104 ( with the results depending on the sample preparation and on the measurement tech¬ 
niques used. Computationally, the band gap is well reproduced with the LDA+U approach (see 
Table [ 2 ]) , however the results depend strongly on the value of U while the HSE06 hybrid functional 
overestimates the gap by approximately 1 eV 60 . 

3.4 CuO surfaces 

The potential use of CuO nanostructures for catalysis, sensing or superhydrophobicity and many 
other application^ 06 -! requires a good knowledge of the surfaces involved in these processes. Only one 
experimental study has been performed to date on CuO surfaces. In this work, it was found, using 
LEED, that the CuO(lOO) surface in UHV conditions does not present a reconstruction 107 . 

DFT studies (using a GGA XC functional with the Hubbard U correction combined, in some 
cases, with ab initio thermodynamics) of the stoichiometric and of some defective surfaces have been 
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Ref. 

lattice (a, b, c ) [A] 

H°\ 

m B [+b] 

Band Gap 
[eV] 

XC 

Penghu 

4.05, 4.06, 5.06 

90.02 

0.0 

0.0 

LSDA 

PengHU 

4.56, 3.27, 4.96 

100.2 

0.63 

1.32 

LSDA+U (U=7.5, J=0) 

Anisimov® 

- 

- 

0.66 

1.9 

LSDA+U (U=7.5, J=0.98) 

DebbichiHU 

4.548, 3.305, 4.903 

99.652 

- 

- 

LSDA+U (U=7.5, J=0.98) 

Ekuma 1 ^ 

4.68, 3.42, 5.13 

90 

0.68 

1.25 

DFT+U (U=7.14) 

Heinemann^h! 

4.588, 3.354, 5.035 

99.39 

0.66 

1.39 

LDA+U (U=7.5, J=0.98) 

^Wu® 

4.55, 3.34, 4.99 

99.507 

0.6 

1.0 

LSDA+U (U=7.5, J=0.98) 

1 RP) 

— 

— 

0.63 

1.1 

GGA+S+U (U=7.5, J=0.98) 

Jiang®°l 

4.68, 3.42, 5.13 

99.54 

0.80 

- 

GGA+U (U=4.5) 

Nolan^ u fl 

4.395, 3.846, 5.176 

- 

0.53-0.7 

0.17-2.11 

GGA+S+U (U=3-9, J=0) 

Svane 102 

- 

- 

0.65 

1.43 

SIC-LSDA 

SzoteklRSI 

- 

- 

0.64 

1.0 

SIC-LSDA 

Heinemann^ 

4.513, 3.612, 5.141 

97.06 

0.54 

2.74 

HSE06 

Lany^ 

- 

- 

- 

1.19 

GW+Vrf 

Exp. 

4.684, 3.423, 5.129 1 ® 

99.54^ 

1 0.68 


0 


Table 2: Calculations of the CuO bulk structure (lattice parameters a, b , c in A, angle /3 (between the 
a and c axes) in degrees, magnetic moment ms in Bohr magnetons and bang gap in eV. The S in the 
functional name indicates spin-polarised calculations. All calculations featuring the Hubbard U (+U 
in the functional name) have been done with the so-called ‘Dudarev approach’ 105 , which requires two 
parameters, U and J. U and J are in eV. 

performed 02 108 109 . CuO(lll) is the most stable surface also compared to defective CuO structures 
with surface and subsurface vacancies 1 ®, with surface energy 7 = 0.74 J/m 2 (although much higher 
than the lowest energy CU 2 O surface^ 7 ). This is true for all oxygen partial pressures, except for a 
narrow range at very low O 2 pressures, where the Cu-terminated CuO(110) surface is favourable. The 
energy order of the stoichiometric surfaces is (111) < (111) < (Oil) < (101) < (110) < (010) < (100) 
(shown in Fig. [3]) 92 . The highest surface energies (2.28 J/m 2 for CuO(lOO)) are associated to large 
fields between anionic and cationic layers, such as for the (010) and (100) surfaces. Among the non- 
stoichiometric surfaces, the O- and Cu-terminated CuO(llO) and CuO(lOO) are more stable than their 
stoichiometric counterpart for high (O-terminated) and low (Cu-terminated) oxygen pressures 92 . 

While important information on stoichiometric and defective CuO surfaces has been obtained using 
computational methods, experimental data are lacking and there is room for work to be performed in 
order to either confirm these predictions or to suggest more possible structures. 

4 Initial stages of oxidation: oxygen adsorption on clean copper 
surfaces 

The growth of an oxide can occur when a metal surface comes into contact with an oxygen-rich 
environment. The reaction sequence leading to oxidation of a clean metal surface is generally accepted 
to be oxygen chemisorption, nucleation and growth of the surface oxide, and bulk oxide growth. 

It is well known that clean copper surfaces in vacuum do not reconstruct 112 . However, exposure 
to oxygen pressure which is low enough not to trigger the formation of oxide-like structures induces 
reconstructions. We will discuss the details of these reconstructions below, and whether they are 
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Figure 3: Side views of CuO slabs showing the structure of stoichiometric cupric oxide CuO surfaces. 
They are shown in order of increasing surface energy (left to right, top to bottom), as calculated using 
ab initio calculations, using the GGA+U method 92 . Copper is shown in brown, oxygen in red. 



Cu(111) 




Figure 4: Overview of the structures formed by oxygen on different copper surfaces, at low oxygen 
exposures, a, d, f) Clean surfaces of Cu(100), Cu(110) and Cu(lll). T, H, B, SB and LB indicate 
respectively the top, hollow, bridge, short bridge and long bridge sites for adsorption, b) c(2 x 2) 
overlayer structure seen on Cu(100) at coverages < 0.3 ML. c) (2\/2 x \/2)i?45° nrissing-row recon¬ 
struction seen on Cu(100) at coverages > 0.5 ML. e) (2 x 1) added-row reconstruction observed on the 
Cu(110) surface, g) Cu20(lll)-like reconstruction of the Cu(lll) surface having hexagonal geometry 
as shown by the dotted line. Two layers of subsurface oxygens are shown. Distorted variations of this 
reconstruction has been observed experimentally 110 111 . 
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important to the initial stages of oxidation is still openly debated. It has been speculated that a 
reconstructed, O-saturated layer must form before the onset of oxidation, since dwell times, i.e. the 
lapse of time between the beginning of oxygen deposition and observation of oxide formation, of up 
to 30 minutes have been observed® 3 115 . Moreover, evidence of the existence of the reconstructed 
copper surfaces and of subsurface growth of the oxide before the onset of island formation, has been 
produced using STM on Cu(100) 116 and Cu(lll) 11 '. DFT calculations have shown that, on Cu(100), 
subsurface oxide-like structures are more easily produced when the surface already has pre-adsorbed 
oxygen atoms, thus suggesting that the reconstruction facilitates oxide formation 118 . However, it 
has also been shown that upon exposure to oxygen at higher pressures 119 copper oxide formation at 
step sites occurs without prior surface reconstruction thus opening the debate about whether direct 
formation of oxide islands can occur without the formation of a O/Cu overlayer. 

In the following section, the very initial stage of the oxidation process, i.e. the chemisorption of 
oxygen onto copper surfaces is reviewed. A summary of the overlayer structures formed on Cu after 
O-dosing is shown in Fig. [4j Cu(100) presents two main reconstructions, one with c(2 x 2) symmetry 
(Fig. [4}}) and a missing-row reconstruction (Fig. [ 4 J 2 ) , which are discussed in depth in Sec. 4.1 The 
main overlayer structure for Cu(110) is the added-row reconstruction shown in Fig. [4j3 and discussed 
in Sec. |4.2[ ), while the 0/Cu(lll) system, presented in Sec. 4.3 shows a distorted hexagonal structure 
resembling the Cu20(lll) surface (Fig. |4^). 

Oxygen adsorption on clean copper has been extensively studied experimentally and computation¬ 
ally, and it has been the subject of a number of reviews over the years, e.g. see Refs. 120 121 . Here we 
focus only on the studies which are relevant in order to understand the onset of copper oxidation on 
the Cu(100), (110) and (111) surfaces and on the most recent developments in the field. 


4.1 Cu(100) 

Oxygen adsorption on Cu(100) has been widely studied® 3 * 121 * 122 -!, and an interesting variety of struc¬ 
tures is formed depending on temperature and coverage, as shown in the phase diagram in Fig. [5^,. 

At low temperatures (up to 100 K) and low coverage (~ 0.1 monolayers, ML, defined as one 
adsorbed oxygen atom for every surface copper atom) experimental and computational evidence 
has shown that incident oxygen molecules dissociate with the oxygen atoms adsorbing at the hol¬ 
low sit d 123 H lS0 These dissociated oxygen atoms stabilize chemisorption of further incoming oxygen 
molecules at higher coverages®®®]. 

Below 473 K, two overlayer structures form: a c(2 x 2) phaseP 6 133 " 136 at a coverage of ~ 0.3 ML 
and a (2\/2 x a/ 2)i?45° missing-row (MR) reconstruction 137 ®® at ~ 0.5 ML. In the c(2 x 2) phase 
(shown in Fig. [ 4 ) 3 ) the oxygen atoms occupy four-fold hollow sites on Cu(100) and form nanometre-sized 
c(2 x 2) domains separated by oxygen-deficient zig-zag shaped boundaries. The MR reconstruction 
(Fig. \$p. [§)) can be viewed as a c(2 x 2) structure with each fourth [100] row of Cu atoms missing. 
Between 0.3 and 0.5 ML, Cu atoms are ejected from the c(2 x 2) domains 118 147 and MR islands start 
forming on terraces, until, at 0.5 ML coverage, a network of missing-row reconstruction islands covers 
the whole surface® 3 133: . 

At high temperatures (473 < T < 1000 K) a c(2 x 2)-like state, with 25% vacancies in the top Cu 
layer forms at 0.5 ML coverage, instead of the MR reconstruction 148 . 

No further oxygen adsorption occurs above 0.5 ML for experiments at very low pressures ( <~ 10 ' 
Torr). However, at higher pressures, further exposure of the surface to O 2 leads to growth of oxide 
in the subsurface regions, while the surface still exhibits the MR reconstruction® 6 * 141 *, validating the 
hypothesis that oxygen-induced surface reconstruction is indeed the first step of oxide growth. 

DFT studies combined with ab initio thermodynamics* 126 * 149 * confirmed that the two experimentally 
observed reconstructions at lower temperatures are the most stable during the early stages of Cu(100) 
oxidation prior to the onset of bulk oxidation. 
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T [K] 

Figure 5: a) Summary of the structures forming on Cu(100) after oxygen exposure. The dashed line 
indicates the lowest temperature at which experiments have been performed (30 K) for this system. The 
areas in white represent domains where no ordered oxygen overlayers form. For low temperatures, at 
coverages of ~ 0.3 ML the c(2x2) structure forms and at ~ 0.5 ML the missing row (MR) reconstruction 
is seen. At intermediate coverages both structures exist. At high temperatures a ‘disordered’ c (2 x 2 ) 
structure occurs. At coverages over 0.5 ML subsurface oxide starts to grow, b) STM image of Cu(100) 
after O 2 exposure (P = 3.7 x ICE 2 mbar, T = 373 K). A (2a/2 x y/2)RAh° island and a zigzag phase 
boundary (bright stripe) are visible. The missing rows of Cu run along the (001) directions and are 
imaged as depressions. Taken from Ref. 116 . 


The transition between the two low-temperature reconstructions has been tentatively explained 
in terms of stress relief, electrostatics and orbital hybridization. Compressive surface stress has been 
shown, both by means of experiments and DFT calculations 150 151 to increase with oxygen adsorption, 
and to be higher in the c(2 x 2) reconstruction. Therefore the MR reconstruction is more stable at 
higher oxygen coverages. Electrostatically, the driving mechanism for oxygen overlayer formation has 
been related to long-range Coulomb interactiorf^2HIi3 and the small size of the c (2 x 2 ) domains to 
repulsion between O adatoms and Cu and O adatoms. At coverages where the O atoms have nowhere 
to form distinct c(2 x 2) domains the phase transition occurs to minimize high 0-0 repulsion. This is 
in disagreement with the findings of Merrick et al . 1 ^, who argue that the stability of the MR system 
is determined by orbital hybridization of neighbouring O-Cu which lowers the energy of the system, 
rather than by long-range interactions. 

The transition between the ordered MR to the ‘disordered’ c(2 x 2) reconstruction has been ex¬ 
plained by means of DFT in terms of diffusion of Cu surface vacancies from an ordered array in the 
MR system to random positions as the temperature increased®. 

An intermediate added row structure was predicted by Kangas et a /.' 159 which has energies compara¬ 
ble to the MR reconstruction. This reconstruction has however not yet been observed experimentally, 
possibly because the right conditions of temperature and pressure for this structure have not been 
used-- 16 l Ab initio thermodynamics calculations could approximately establish at which experimental 
conditions this reconstruction is expected and could thus inform further experimental work. 

It is important to note that simulations have shown that the MR structure is a necessary step 
towards the formation of a Cu 20 -like structure 118160 . Indeed, Cu 20 -like structures were found to 
form on the MR reconstruction upon O adsorption which were not found on the non-reconstructed 
surface, thus confirming the experimental results of Lahtonen et al . 116 and fitting well with the work of 
Zhou et al . 161 , who finds the presence of a two-layer thick oxide before the formation of oxide islands 
of the Cu(100) surface (see Sec. 5.1). 


11 





















In summary, three reconstructions can form on Cu(100) upon oxygen adsorption, as shown in 
Fig-! Of these, the MR reconstruction is the most stable at room temperature, for coverages ~ 0.5 
ML. Experimental and computational evidence has shown that the formation of this reconstruction is 
the first step towards the formation of the bulk oxide. 

4.2 Cu(110) 

Molecular oxygen dissociates when deposited on Cu(110) at temperatures above 45 K 162 . Upon dis¬ 
sociation of the O 2 molecules a number of overlayer structures are observed, including an added-row 
(2 x 1) structure (Fig. [ 4 J 3 ) which is formed at an oxygen coverage of 0.5 ML 121 145 163 165 . Another 
type of surface reconstruction, c(6 x 2 ), has also been reported at high coverage ( 2/3 MLjEsrami 

or at lower coverage but higher temperature 173 . 

The (2x1) phase forms via the creation of Cu-O-Cu-O chains along the [001] direction of the 
substrate, which eventually become the ‘added rows’ on top of the clean substrate at a coverage of 
0.5 ML. The chains start forming above 70 K but do not fully organize until ~ 200 K 162 . They are 
formed from mobile chemisorbed O atoms and Cu adatoms which leave from step edges and diffuse 
across the terraces 138 162 174 H 178 . The experimental barrier calculated for the formation of these strings, 
0.22T0.01 eV^ 179 , is close to the DFT-calculated barriers for Cu (0.25 eV) and O (0.15 eV) diffusion! 180 ! 

At oxygen coverages between 0.05 and 0.45 ML these Cu-O-Cu-O chains self-organize in a periodic 
array (called a supergrating) with a spacing varying between 60 and 140 A 138 181 184 . The dependence 
of the period of the supergrating as a function of oxygen coverage has been explained in terms of 
electrostatic 185 and elastic interactions 186 191 . In the first case, the period has been related to the 
difference in work function between the clean and reconstructed sections of the surface, in the second 
case to stress relief from the mismatch between the preferred period of the Cu-O-Cu-O chains and the 
period of the Cu(110) substrate in the [001] direction. It could indeed be possible for both mechanisms 
to be at play. 

The stability of the (2 x 1) reconstruction at 0.5 ML has been established computationally by 
means of semiempirical and DFT calculations, where it has been found to be energetically favourable 
with respect to the unreconstructed surface 179 180 192 19 ® and to an alternative added row (4 x 1) ge¬ 
ometry 190 . DFT simulations with ab initio thermodynamics 126 1,2 191 found that the Cu-0 added row 
reconstruction is favoured at low oxygen exposures, whereas at higher oxygen exposures, a transition to 
the c(6 x 2) structure is predicted to occur. Thus, experiment and theory are in qualitative agreement, 
although the absolute transition pressures vary enormously (by 10 orders of magnitude). Since it is 
a major challenge for current DFT XC functionals to accurately predict adsorption energies and the 
underlying value of the adsorption energy has a huge impact on subsequence pressure estimates, such 
a quantitive discrepancy is not uncommon • ■ • . The barrier to transition between the two structures 
has been calculated to be fairly high at 1.41 eV, which seems to explain why this phase is observed 
only at high temperatures. 

A large amount of work has been performed in order to understand the surface electronic states of 
the added row reconstruction, using both experimental methods such as photoemission spectroscopy 195 198 
and theory! 198 200 . The character of the O-Cu bonding is found to be predominantly ionic, and the 
surface 0(2 p) orbitals hybridize strongly with the Cu(3d) states, forming bonding and antibonding 
linear combinations, with the antibonding bands not having been identified unambiguously yet. 

There is no evidence so far on how the observed reconstruction relates to the initial formation of 
the oxide, if at all. DFT calculations of subsurface oxygen added beneath both reconstructed surfaces 
have found that when an O coverage of 1 ML is reached, subsurface oxide formation in the tetrahedral 
interstitial sites is predicted to occur 201 . The presence of oxygen in the tetrahedral site has been linked 
to oxide formation (since the O in CU 2 O resides in the tetrahedral sites), and therefore this is a possible 
mechanism for the initial formation of the oxide, which should however be confirmed by experimental 
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or further computational work. 


4.3 Cu(lll) 

The clean (111) surface has the lowest surface energy for copper and it is less reactive compared to the 
other low-index Cu surfaces. It is therefore less studied for oxygen adsorption ^ 10 111 114 129 130 202 "2 10 '. 
No ordered structures are observed experimentally for low oxygen exposure 114 201 205 211 . Indeed, also 
DFT studies have found that oxygen adsorbs preferentially at the threefold hollow site for coverages 
up to 0.75 ML, without forming periodic overlayer structures 49 212 214 . 

At higher coverage, the adsorbed oxygen and copper adatoms that are ejected from step edges 
and terrace d 117 I 209 I 2 IS] start forming overlayer structures, before the onset of epitaxial growth (see 
Sec. 5.3). Two main classes of reconstructions have been proposed for the 0/Cu(lll) surface. The 


first, seen in LEED studies, involves a Cu(100)-like overlayer, incommensurate with the underlying 
unreconstructed Cu(lll) surface, with the oxygen atoms occupying the hollow sites 215 216 . The second 
class comprises two long-range ordered structures 11 0 l 207 l, the so-called ‘29’ and ‘44’ superstructures. 
They have, respectively, (Vl37?(46.1 0 x 7R21.8°) and (V73-R5.8° x \f2lR — 10.9°) symmetry and very 
large surface unit cells, 29 and 44 times larger than the lxl cell of clean Cu(lll). They exhibit an 
honeycomb pattern formed by distorted honeycomb units of the ‘ideal’ Cri 20 (lll) overlayer shown in 
Fig- The oxide-like overlayer structure has been confirmed by STM 2 ®, X-ray absorption studies ^ 1 ^ 
and DFT and ab initio molecular dynamics 49 ' 8 . This differs from the behaviour of other O/transition 
metal systems such as 0/Ag(lll) 21 ' 218 and 0/Pd(lll) 219 , where chemisorbed oxygen adlayers form, 
before the formation of a surface oxide. Other systems, such as 0/Ru(0001 ) 220 and 0/Rh(lll ) 221 2 ®, 
instead never form a surface oxide layer before the onset of bulk oxidation. 

Therefore, both experimental and computational studies point to the formation of hexagonal or 
quasi-hexagonal structures on the Cu(lll) surface upon oxygen adsorption, structures which can be 
viewed as the initial layer of a Cu 20 (lll) film and which can potentially act as a template for the 
growth of further Cu 20 (lll) layers. 


5 Oxide film growth 

When copper surfaces are exposed to a continuous flow of oxygen for a long time oxidation is expected 
to occur. It is known in general, from well-defined surface science studies of metals, such as those 
discussed in Sec.[4j that the O/metal structure which forms at the onset of oxidation can be distinct 
from the bulk oxide surface structure 14 223 224 . The growth of oxides in copper has been studied since 
the early 1920s and early work has been extensively reviewed (see e.g. Ref. 11 ). In this section, we aim 
to present an overview of the main topics in the field of copper oxide growth mainly focusing on recent 
studies and on the present state of experiments and theory. 

In order to understand the microscopic details of the onset of oxidation, experiments have been 
performed in controlled laboratory conditions, where the orientation of the surface, temperature and 
oxygen pressures can be tuned to the required values. The formation of the oxide has been directly 
observed by means of TEM, and their structure analysed by means of LEED and ellipsometry. At 
the oxygen pressur e, te mperature and exposure times used in these studies, only cuprous oxide is 
expected to form 39 225 226 . We review studies showing that the growth of cuprous oxide on low-index 
copper surfaces is epitaxial with the substrat e) 1 4 l 227 H 25ffl through nucleation and coalescence of nano¬ 
islands 113 114 202 203 227 —®232] ^ ag sc h em atically represented in Fig. [ 6 ]). However, the mechanism of 
formation of the islands and the resulting shapes and growth rates are strongly dependent on the 
Cu substrate, as shown in Sec. 5.1 for Cu(100), Sec. |5.2| for Cu(110) and Sec. 5.3 for Cu(lll). Key 


experimental findings are summarised in Tables iiH The kinetic models for the nucleation and 
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Figure 6: Schematic illustration of the stages of growth of a copper oxide film on copper, a) Gas-phase 
oxygen molecules and clean copper surface, b) Upon adsorption of the oxygen, the copper surface 
reconstructs, c) Nucleation of oxide islands upon diffusion of oxygen on the reconstructed surface, 
until the island saturation density ( N sa t ) is reached, d) Growth of the oxide islands. As the islands 
grow bigger direct oxygen impingement on the islands starts playing a more important role, e) After 
the islands coalesce oxide growth proceeds through interfacial diffusion of oxygen. 


coalescence of the nano-islands and for the growth of the oxide thin film are reviewed in Sec. 5.4 and 
Sec. 15.51 

Growth of the oxide in ambient conditions, when the copper surface is exposed to humid air, has 
also been extensively studied. The body of work addressing this issue is reviewed in Sec. 5.5 and 
summarised in Table [6l 

Finally, the existence of preferential oxidation sites is discussed (Sec. 5.7) and studies on the initial 
growth of CuO are reviewed (Sec. 5.8). 


5.1 Oxide nano-islands: Cu(100) 

As for the formation of O/Cu overlayers, Cu(100) is the most extensively studied surface, both ex¬ 
perimentally and computationally. Copper oxidation on this surface proceeds, for low oxygen partial 
pressures (P ~ 10 -4 Tori'), through the formation of islands which are epitaxial with the surface, i.e. 
Cu 2 O( 100 ) || Cu(100/ 113 * 115 * 231 l ^^^ and have 6x7 lattice misfit configuratkm)225I232 These islands 
grow and coalesce with further oxygen deposition. The shape of the islands depends on temperature 238 
in a rather interesting manner as shown in Fig. [7J Below 350 °C only triangular islands are observed, 
whereas between 400 — 550 °C the shape changes to round or square. At around 600 °C the islands start 
to grow until, at a critical size of ~ 110 nm, when they switch to a quasi-one-dimensional elongated 
rod shape. Between 650 — 800 °C pyramid-shaped islands are observed and between 800 — 1000 °C 
hollow pyramids form 239 . The big effect of temperature on the morphology of the islands could be 
either due to the dependence of copper and oxygen diffusion on temperature or to changes in interfacial 
strains and in the mechanical properties. Indeed, the transition between the round/square to rod shape 
shown in Fig. 0h c has been explained as a competition between the surface energy and the strain of 
the islands due to the mismatch between the clean metal and the oxide lattices, following the so-called 
Tersoff-Tromp 240 energy model. 

TEM, XPS and AES studies looking at the cross section of an oxidising Cu(100) single crystal under 
low oxygen partial pressure 115 161 — 242 have revealed that the oxide islands form on top of an oxide 
wetting layer. The wetting layer itself has a (\/2 x 2\/2)i?45 0 missing row reconstruction and forms 
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Figure 7: Oxide islands observed on the Cu(100) surface using in situ TEM imaging, for a set of 
temperatures between 350 — 1000 °C. The oxidation partial pressure is P = 5 x 10~ 4 Tor. A large 
variety of island structures are observed, according to the oxidation temperature. Taken from Ref. 238 . 


4.1) who observed 
~ 10 ^ Torr. 


from the Cu atoms ejected by the missing row reconstruction of the substrate (described in Sec. 4.1) 
and the further oxygen deposited on the substrate!® . Island growth also occurs beneath the surface, 
in good agreement with the already reviewed results by Lahtonen et al. (Sec. 
subsurface oxide growth for oxygen coverages above 0.5 ML and pressures above 

In general, when a metal is capable of forming uniform subsurface oxides (as in the case of 
Ag(110) 243 and Ru(0001)^®), oxide growth proceeds uniformly, rather than via island formation. 
The reason for the island formation on top of the wetting layer is found in stress mismatch between 
the Cu(100) substrate and the Cri20(100) film 161 . 

The transition between the O/Cu(100) system described in Sec. 4.1 and the oxide islands has 
been studied to some extent with DFT, and although the full transition has not been modelled yet, 
insight has been gained into the mechanisms at play. On the MR reconstructed surface, sub-surface 
adsorption becomes favourable for O coverage above 1.0 ML, with low barriers for the transport of 
oxygen atoms below the surface^- 24 ' 5247 . Sub-surface oxygen atoms below the missing-row reconstruction 
adsorb in the tetrahedral interstitial sites, and thus form an oxide-like structure, unlike the case of 
sub-surface oxygen atoms below non-reconstructed surfaces, which adsorbs in the octahedral site 16 ®. 
However, the limiting factor for the oxidation of Cu(100) is the dissociation and on-surface diffusion 
of the oxygen molecule. The dissociation of the oxygen molecule, although almost barrierless on the 
clean Cu(100) 12 ^, is blocked by the on-surface oxygen on the reconstructed surfacd 248 , and requires 
the diffusion of the oxygen molecule towards either a high-Cu concentration ared 1 ' 53 or vacancies and 
surface defects 249 250 where the dissociation barrier is lower. Diffusion of O and Cu atoms on the 
MR reconstructed surface is slow, having barriers of 1.4 eV for oxygen and 2.0 eV for Cu 153 . The 
high barriers and slow diffusion make these processes difficult to simulate. Indeed, Devine et al. 251 
successfully simulated oxygen molecule dissociation on the clean Cu(100) surface with bond-order 
potentials, however they did not see any dissociation event on the missing-row reconstructed surface 
within their (short) 10 ps of molecular dynamics. This highlights the difficulty of reproducing these 
complex phenomena which include diffusion of oxygen molecules, dissociation, diffusion of oxygen and 
copper atoms on and through the surface using standard computational approaches. To this end, the 
combination of an accurate reactive potential with methods which allow for rare events to be probed, 
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such as 


metadynamics 252 are possibly a way to go for this type of system. 


Reference 

Sample 

thickness 

Surface Preparation 

Exp. Condi¬ 

tions 

Technique 

Result 

Brockwayl^hJ 

(1972) 

90 nm 

Annealing in H2 at T=630 
°C 

P=10 -ii Torr, 

T=525 °C 

TEM 

Epitaxial growth of ox¬ 
ide islands. 

Heinemannhh* 

(1975) 

1 80 nm 

Annealing, argon ion sputter 
etching 

P=5 x 10“ a Torr, 

T=425 °C 

TEM 

Epitaxial growth of 
square and hexagonal 
oxide islands. 

Stefanov^ - 

(1988) 


Ion bombardment and heat¬ 
ing under UHV 

10 to 2 x 10° L, 
T=—130—180 °C 

HREELS, 

XPS 

CU2O formation at oxy¬ 
gen exposures between 
10 5 - 2 x 10 6 L 

Yang‘S 

(1997) 

40 nm 

Annealing in methanol 
vapour at 350 °C 

1.5 x 10“ 5 Torr 

TEM 

Rate of growth of the ox¬ 
ide islands oc t 1 ' 3 . 

Yang^ 

(1998) 

lOOnm 

Annealing in UHV/CH4O 
vapour at 350 °C 

1 x 10~ & — 1 x 10~ 4 

Torr 

TEM 

O monolayer forms be¬ 
fore growth of the oxide 
islands. 

YangP^ 

(1998) 

60 nm 

Annealing in methanol 
vapour at 350 °C, 5 x 10~” 5 
Torr 

5x 10 -5 —760 Tor 
at 60 - 600 °C 

TEM 

Epitaxial island forma¬ 
tion after surface recon¬ 
struction. 

Yang^ 

(1998) 

100 nm 

Annealing in CH4O vapour 
at 350 °C, 5 x 10 -5 Torr 

5x10~ 4 Tor 290- 
435 °C 

TEM 

Nucleation of islands 
promoted by O diffusion. 

Yang‘S 

(1999) 

100 nm 

Annealing in CH4O vapour 
at 350 °C, 5 x 10 -5 Torr 

Ix5“°-5xl0' 4 
Tor 70 - 600 °C 

TEM 

Preferential nucleation 
site at the edge of holes. 

Yang*^l 

(2001) 

60 - 100 

nm 

Annealing in CH4O vapour 
at 350 °C, 5 x 10 -5 Torr 

O2/H2O vapour 
at 5 x 10 -4 Tor at 
350 °C 

TEM 

Initial surface recon¬ 
struction prior to island 
growth. 

Yan^ 

(2002) 

60 nm 

Annealing in CH4O vapour 
at 350 °C, 5 x 10 -5 Torr 

5 x 10' 5 - 760 

Torr at 60 — 600 
°C 

TEM 

Good agreement of 
kinetic data with the 
JMAK model. 

Zhou 1 ^ 

(2003) 

70 nm 

Annealing in CH4O vapour 
at 350 °C, 5 x 10 -5 

5 x 10“ 4 Torr at 
150 - 1000 °C 

TEM 

Temperature-dependent 
shape of oxide nano¬ 
islands. 

Eastman^" &: 

(2005) 

110 - 200 

nm 

Annealing in Ar-2%H2 at 
850 °C, 5 x 10 -5 Torr 

5 x 10 -4 Torr at 
350 - 780 °C 

X-ray 

scattering 

Epitaxial nano-island 

formation. 

“Zhou® 

(2005) 

70-80 nm 

Annealing in CH4O vapour 
at 350 °C/vacuum 800° C 

5 x 10~ 4 Torr, 350 
°C 

TEM 

Island nucleation rate 
faster on (111) than 
(110) or (100) 

Zhou 1 ^ 

(2005) 

70-80 nm 

Annealing in vacuum at 550 
°C 

5 x 10“ 5 Torr, 

350 - 900 °C 

TEM 

Temperature-dependent 
shape and oxidation rate 
of the oxide islands. 

Zhou 1 ^ 

(2005) 

70-80 nm 

Annealing in CH4O vapour 
at 350 °C/vacuum 800 °C 

5 x 10~ 4 Torr, 350 
°C 

TEM 

Island nucleation rate 
faster on (111) than 
(110) or (100) 

Lampimaki 1 ^ 

il 

Ar + bombardment, anneal¬ 
ing at 700 °C 

2.8 x 10~ 2 - 160 
Torr, 30- 100 °C 

XPS, 

XAS, 

STM 

Island formation on top 
of the missing-row recon¬ 
struction for O/Cu(100) 

Lahtonenhl^l 

(2008) 

0.5 mm 

Ar + bombardment, anneal¬ 
ing at 430 °C 

6 x 10“' Torr and 
2.8 x 10~ 2 Torr, 

T = 100 °C 

STM 

Surface reconstruc¬ 

tion and oxide island 
formation at high O2 
exposure. 

Zhou 1 ^ 

(2009) 

70 nm 

Annealing in Ar/H2 at 700 
°C 

5 x 10“ a Tor, 700 
°C 

TEM 

CU2O islands, 200/500 
nm side. Epitaxial 

growth with inclined 
CU2O/CU edges. 

Zhou 1 ^ 

(2012) 

50 nm 

Annealing in H2 at 600 °C 

5 x 10“° Tor, 350 
°C 

TEM 

Step-edge induced oxide 
growth. 
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Zhou 6 ^ 

50 nm 

Annealing in H 2 at 600 °C 

1 x 10 _il Tor, 550 

TEM 

CU 2 O islands grow on an 

(2013) 



°C 


oxide wetting layer nu¬ 
cleated on surface steps. 


Table 3: Experimental results of copper oxide formation in ultra-high 
vacuum on Cu(100). Experimental conditions and techniques have been 
listed, together with a summary of the main result of each study. 


5.2 Oxide nano-islands: Cu(110) 


-4 


The islands formed on Cu(110) single crystals after being exposed to low pressure oxygen (P ~ 10 
Torr) are also found to be epitaxial with the substrate 182 . The island morphology and the time required 
to reach saturation density depend on temperature, with higher temperatures requiring much shorter 
time to saturation. Between 450 and 650 °C the lateral size of the islands is almost constant (200 — 250 
nm), however the islands were found to thicken considerably (24 — 40 nm) beneath the Cu surface, 
showing that the rise in temperature greatly enhances the interfacial diffusion of oxygen 25 '. At 700 
°C, the (110) surface of clean copper roughens with a step height of ~ 10 nm 258 . Exposure to oxygen 
of this rough surface shows the formation of a higher density of oxide islands with a fast nucleation 
rate but slower lateral growth than for smooth surfaces. This comparison between smooth and rough 
surfaces supports the idea that the kinetics of nucleation and growth of islands is dependent on the 


surface diffusion of oxygen atoms on the surface (see Sec. 5.4) 


Reference 

Sample 

thickness 

Surface Preparation 

Exp. Condi¬ 

tions 

Technique 

Result 

^ZhorPU 

( 2002 ) 

70 nm 

Annealing in CH 4 O vapour 
at 350 °C, 5 x 10 - 5 Tor 

8/1 x 10 -4 Tor at 
600 °C 

TEM 

Rod-shaped 20 nm-thick 
islands. 

“ZhorP 3 

(2003) 

70 nm 

Annealing in CH 4 O vapour 
at 350 °C 

5 x 10 “ 4 Tor, 300 
°C and 450 °C 

TEM 

Faster initial oxidation 
rate than Cu(100) 

“ZhorP 3 

(2004) 

70 nm 

Annealing in CH 4 O vapour 
at 350 °C 

5 x 10 “ 4 Tor, 350 
°C and 750 °C 

TEM 

Higher density and 
slower lateral growth 
rate of islands on 
rougher Cu(110) sur¬ 
faces. 

^ZhorP 3 

(2004) 

70 nm 

Annealing in CH 4 O vapour 
at 350 °C 

5 x 10 “ 4 Tor, 750 

°C 

TEM 

Faster oxidation at 

higher temperatures. 

“ZhorPSI 

(2005) 

70-80 nm 

Annealing in CH 4 O vapour 
at 350 °C/vacuum 800° C 

5 x 10 “ 4 Tor, 350 

°C 

TEM 

Island nucleation rate 
faster on ( 111 ) than 
( 110 ) or ( 100 ) 

~ZhcnpSl 

(2005) 

70-80 nm 

Annealing in vacuum at 550 
°C 

5 x 10 -5 Tor, 

350 - 900 °C 

TEM 

Temperature-dependent 
shape and oxidation rate 
of the oxide islands. 

“ZhorP 3 

(2005) 

70-80 nm 

Annealing in CH 4 O vapour 
at 350 °C/vacuum 800 °C 

5 x 10 “ 4 Tor, 350 
°C 

TEM 

Island nucleation rate 
faster on ( 111 ) than 
( 110 ) or ( 100 ) 


Table 4: Experimental results of copper oxide formation in ultra-high 
vacuum on Cu(110). Experimental conditions and techniques have been 
listed, together with a summary of the main result of each study. 


5.3 Oxide nano-islands: Cu(lll) 

At room temperature the oxidation of Cu(l ll) also proceeds through the formation of epitaxial oxide 
islands with a 5 x 6 lattice misfit 11 ' 209 260 which coalesce with continuous oxygen exposure. Three 
processes of oxide formation are observed in this regime 117209 (Fig. [8]): growth from step edges, in¬ 
terrace growth (~ 1.2 A below the terrace surface) from vacancy islands and growth of on-terrace 
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oxide. The formation of randomly placed islands at low temperature agrees well with the existence of 


a disordered underlying structure observed for low-coverage oxygen adsorption and analysed in Sec. 4.3 


At higher temperatures 209 260 fast formation of the oxide and fast lateral growth leads to the formation 
of a two-dimensional oxide structure. This uniform growth is templated by the distorted-hexagonal 
surface reconstruction shown in Sec. ER 

Zhou et al. observed an interesting phenomenon at intermediate temperatures: the islands nucleate 
close to existing islands, anisotropically elongating along the [110] direction in a percolating manner, 
as shown in Fig. [9] 



Figure 8: Structures observed with STM by Matsumoto et al.- 209 when oxidising a Cu(lll) surface with 
step edges (visible in panel a). In-terrace and on-terrace oxides are visible in panel b, and atomically 
resolved in panel c. In-terrace oxide is darker and on-terrace oxide lighter than the clean copper 
surface. Panel d shows the oxide which grows at the edge of the terrace. The area show in panel a,b is 
2000 x 2000 A 2 , in panel c,d is 200 x 200 A 2 . The square in panel a indicates part of the area shown 
in b. Squares and the arrow in b indicate the areas shown in c and d. 

This ‘discontinuous-branched’ shape has been investigated using kinetic Monte Carlo techniques^H, 
and it appears to be related to restricted diffusion of oxygen on the surface, which might then be related 
to the surface being reconstructed to the ‘29’ or ‘44’ structures, which are fairly corrugated (they extend 
up to 3.1 A over the clean Cu(lll) surface). 


Reference 

Sample 

thickness 

Surface Preparation 

Exp. Condi¬ 

tions 

Technique 

Result 

Lawles#^ 1 

(1956) 


Electropolishing 

P = 0.8 - 760 
Torr, T = 170 - 
450 °C 

XRD 

Epitaxial CU 2 O oxide, 
disordered CuO. 

Goulden^ 3 ” 

(1976) 

250 - 400 
/mi 

Polishing in orthophospho- 
ric acid at 1.4 V 

8 x 10“ 4 Torr, 

T = 250-400 °C 

TEM 

Epitaxial oxide islands. 

(1978) 

— 

Annealing in vacuum at 300 
°C 

10“ 4 -10“' Torr, 

T = 350°C 

TEM 

Epitaxial lamellar 

growth. 
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Dubois 1 ^ 

(1982) 

— 

Annealing in H 2 at 630 °C 

1 x 10 3 Torr, 

T = 525 °C 

EM 

Epitaxial oxide islands 
on defect sites. 

Milne! 2 ^ 

(1984) 

0.05 fim 

Polishing, annealing 

P = 10 -b - 5 x 
10 -4 Torr, T = 
300 °C 

RHEED, 

TEM 

Formation of epitaxial 
oxide islands. 

Rauh® 

(1993) 

50 nm 

Deposition by dc heating. 

0 2 , P = 7.5 x 
10 -4 - 9 x 10 -2 
Tor, T = 105 °C, 
400 min 

Ellipsom. 

Formation of a CU 2 O 
film 

Matsumoto 1 ^ 

( 2001 ) 


Ar + sputtering and vacuum 
annealing at 500 °C 

P = 10 -Y - 10'° 
Torr, RT 

STM, 

LEED 

Growth of oxide from 
step edges. 

Zhou 1 ^ 1 

(2005) 

70-80 nm 

Annealing in CH 4 O vapour 
at 350 °C/vacuum 800° C 

5 x 10 -4 Tor, 350 
°C 

TEM 

Island nucleation rate 
faster on ( 111 ) than 
( 110 ) or ( 100 ) 

-zho^m 

(2005) 

70-80 nm 

Annealing in vacuum at 550 
°C 

5 x 10 -b Tor, 

350 - 900 °C 

TEM 

Temperature-dependent 
shape and oxidation rate 
of the oxide islands. 

Zhou 1 ^ 1 

(2005) 

70-80 nm 

Annealing in CH 4 O vapour 
at 350 °C/vacuum 800 °C 

5 x 10 -4 Tor, 350 
°C 

TEM 

Island nucleation rate 
faster on ( 111 ) than 
( 110 ) or ( 100 ) 

ZhoJSEI 

(2008) 

— 

Annealing in vacuum at 800 
°C 

3 x 10 -4 Tor, 900 
°C 

TEM 

Growth of oxide hollow 
pyramidal islands. 

Zhou 1 ^ 1 

(2009) 

60 nm 

Annealing in CH 4 O vapour 
at 350 °C/vacuum 800 °C 

5 x 10 -4 Tor, 350 
°C 

TEM 

Epitaxial CU 2 O ~ 2.5 
nm-high islands. 

Leoii^ 

( 2012 ) 


Ar + sputtering and 550 °C 
heating 

P = 10~ Y Torr, 

RT, pulse injec¬ 
tion of air 

AES, 

STM 

Oxide growth from step 
edges. 


Table 5: Experimental results of copper oxide formation in ultra-high 
vacuum on Cu(lll). Experimental conditions and techniques have been 
listed, together with a summary of the main result of each study. 


5.4 Nano-island formation kinetics 


A number of models have been proposed to explain the kinetics governing the initial stages of oxide 
growth, from island nucleation to their coalescence. 

The mechanism of nucleation and initial growth of oxide islands on a clean Cu surface has been 
proposed in terms of ‘capture zones’, a well-established idea in non-homogeneous film formation the¬ 
ory 264265 . In this interpretation, schematically shown in Fig. 10 1 , oxygen atoms landing in a capture 
zone, i.e. an area with radius Ld surrounding the island, will likely aggregate with the island and 
contribute to island growth rather than nucleating new islands. Large ‘capture zones’ (i.e. large val¬ 
ues of L(j) are associate with longer path lengths of surface oxygen diffusion. In this framework, the 
saturation number of islands that can be nucleated, N sat , is expressed as" 


N sat = -^(1 - e~ kL ^) 

L d 


(1) 


where k is the initial island nucleation rate and t is time. This model was found to fit well experimental 
data for oxidation of Cu(100) and Cu(110) 234 2 42, as shown in Table [5^4} The smaller number of islands 



N sat [/till -2 ] 

k [/im 2 min x ] 

L d 

(100) 

0.83 

0.17 

1.09 

(110) 

4.34 

1.74 

0.33 


nucleating on Cu(100) is related to a larger oxygen capture area, i.e. longer path lengths for surface 
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Figure 9: Percolating oxide growth on Cu(lll) observed using in situ TEM imaging at 450 °C. The 
oxidation partial pressure is P = 5 x 10 -4 Tor. The layer is formed by preferential nucleation of 
islands along the [110] direction, suggesting that between 350 — 450 °C the substrate transitions from 
a disordered to an ordered structure. Taken from Ref. 256 . 


oxygen diffusion. This is consistent with the structure of the missing-row reconstructed Cu(100) surface, 
which is fairly smooth (it has a corrugation of 0.35 A 138 ), see Sec. 4.1 thus favouring surface diffusion. 


Similarly, the short diffusion path length for Cu(110) is consistent with the added-row reconstruction 
which protrudes ~ 1.5 A over the Cu surfaced 21 * 147 * (see Sec. 4.2). Nucleation is seen to be faster 


on Cu(110) than on Cu(100), see Table 5.4, and even faster on Cu(lll) at 350 o C' 256 } The O-induced 
reconstruction of the Cu(lll) surface has been reported to have a disordered surface overlayer at room 
temperature, with O and Cu atoms at different height^ 10 * 207 ! (see Sec 4.3). This could explain a shorter 
path length for oxygen diffusion on Cu(lll) leading to a fast nucleation of a large number of oxide 
islands. 

After island nucleation has reached saturation point, they start growing until coalescence. Is¬ 
land growth models based on oxygen impingement and surface diffusion (schematically represented in 
Fig. 10 3 ) have been developed 233 266 267 : oxygen surface diffusion initially dominates the oxide growth, 
and later oxygen direct impingement becomes significant when the oxide islands grow larger in size. In 
particular the model proposed by Yang et al . 233 was found to fit well experimental data on Cu(100), 
with island growth proportional to t 1 ' 3 . 

Another well established theory for the formation of thin films is the Avrami (or JMAK) nucle¬ 
ation model 268 " 270 . It presumes isotropic and homogeneous nucleation of the islands and depends 
exponentially on time: 

X{t) = 1 - e~ ktn , (2) 


where X is the oxide thickness. Cu(100) fits this model 242 with k = 1.9 x 10 -4 and n = 2. The value 
of k is much smaller than expected (typically it is k = 7r/3), possibly because of the non-constant 
island nucleation in cuprous oxide (but rather following the relation in Eq. [l]) and the non-linear island 
growth rate (A(t) oc t 1 ' 3 ). Yang et al. modified the JMAK model in order to take into account these 
two factors and found an excellent fit with Cu(100) experimental datd 235 . 

Although the nucleation of islands is faster on Cu(110), the long term (> 60 minutes) oxidation of 
the (100) surface of copper is much faster than on the (110) and (111) I 228 I 23 °I2 >U surfaces. Fast initial 
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Figure 10: a) Schematic representation of the mechanism of island nucleation and growth according to 
the model in Eq. [l] The dashed circles represent the ‘oxygen capture zones’ of radius around the 
oxide islands. Oxygen atoms (in red) diffuse on the surface until they enter a capture zone, when they 
are incorporated in an oxide island (represented as a brown triangle), b) Schematic representation of 
the mechanisms of island growth to coalescence. As the islands increase in size their areas of oxygen 
capture increase (black dashed line). Oxygen surface diffusion (green arrows) is still important, however 
direct impingement of oxygen atom in the islands is also more likely to occur, c) Oxidation of Cu(100), 
Cu(110) and Cu(lll) copper surfaces using the ReaxFF force field at T = 300 K and T = 600 K. 
Cu(lll) is found to oxidise more easily than (110) and (100) and form a thicker oxide film. The 
amorphous nature of the formed oxide is clearly visible. Taken from Ref. 2 ' 1 . 


nucleation and growth of islands on Cu(110) and Cu(lll) leads quickly to a thinner coalesced film, 
which then continues growing through oxygen diffusion through the oxide layer, a much slower process 
than surface diffusion. On the contrary, the slow nucleation of islands on Cu(100) leads to slower 
coalescence to a thicker oxide film with respect to the Cu(110) and Cu(lll) and thus quicker oxide 
growth by means of surface diffusion. 

The case of oxidation on Cu(lll) at higher temperatures (over 550 °C) is the only one where CU 2 O 
two-dimensional thin film growth is observed 117 256 . The O-induced surface reconstruction of Cu(lll) 
at high temperatures, as seen on Sec. 4.3, has a hexagonal morphology similar to the Cu20(lll) 
plane 110 207 which can act as a ‘template’ structure for the layered growth of the oxide. 

A few computational studies have tried to obtain further insight into the kinetics of oxide growth, 
with limited success so far. The correct relative rates on the three low-index Cu surfaces have been 
obtained using molecular dynamics with ReaxFF, a bond order potential 2 ' 1 . However the calculation 
led to the formation of a uniform amorphous thin film (as shown in Fig. [To];), thus not reproducing 
crystalline island formation as seen in experiments, and the role of the surface reconstruction was 
not taken into account. As already mentioned, Devine et a/. 251 , looked to reproduce oxidation on an 
O-reconstructed Cu(100) surface, using molecular dynamics and a bond order potential. However, the 
process was too slow and the simulation too short for it to be modelled. 
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Figure 11: Overview of growth models for copper oxide present in the literature. The graphs show 
the oxide film thickness vs oxidation temperature for a number of experimental studies. The data 
have been divided into two categories, according the composition of the initial copper sample, single 
crystal or polycrystalline. They have also been colour coded according to the oxidation kinetic model 
which has been attributed to them. It can be seen that on the basis of these three parameters only 
(temperature, film thickness and initial sample structure) a single kinetic model cannot be identified. 

5.5 Long-term copper oxidation 

After nucleation, growth and coalescence of copper oxide islands, oxidation proceeds through the 
diffusion of oxygen atoms through the oxide and the character of the oxidation process changes. 

The kinetics and mechanism of copper oxidation, after the exposure of the copper surface to oxy¬ 
gen flow for an extended period (up to 32 hours) have been extensively studied ^ 4 * 272 * 273 !, however no 
consensus on a single rate law describing the oxidation dynamics has been established. 

Clean copper surfaces have been found to oxidize at different rates, with the Cu(100) face reported 
to have the fastest rate of oxidation, at odds with the oxidation rates measured for the formation and 
growth of oxide nano-islands. Indeed, Young 22 ^ et al. and Gwathmey 22 ® et al. found that the order 
of oxidation rate for the low-index surfaces is (100), (111), (HO) with (100) the fastest oxidising facet, 
for a wide range of temperatures. Rhodin 2 ' 4 found instead the order to be (100), (110), (111) with 
(100). This difference in ordering could be due to several factors, for example the use of very different 
experimental analysis techniques. Gwathmey 22 ^ et al. used diffuse white light to analyse the oxidized 
surfaces, observing different interference colours corresponding to different oxide thicknesses, Young 228 
et al. used ellipsometry which relies on a ‘guess’ of a surface model and Rhodin^H calculated the film 
thickness from weight changes in the sample. All three methods present some drawbacks (diffuse white 
light does not provide a value for the thickness and relies on visual checking of the surface colour, 
ellipsometry relies on a guess model and Rhodin’s method depends on an accurate estimation of the 
effective area of the resulting oxide) and it is difficult to tell which of the three is more reliable. The 
difference observed could also be the result of environmental factors such as the presence of impurities or 
different sample preparations, which have been shown to lead to different resulting metal oxides® 73 27 ® 
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which makes oxidizing metals very challenging to study. 

Very large discrepancies are observed in models for oxidation rates, i.e. growth of thickness of 
the oxide as a function of time, and they are summarised in Fig [lTJ A number of theories have been 
proposed®^®®*, mostly based on the Cabrera and Mott 281 theory and they postulate that, under the 
assumption of uniform epitaxial growth, the thickness of a metal oxide increases following an inverse 
logarithmic rate law for very thin films (up to 7.3 nm for Cu) and a cubic law for thicker films (up 
to 1.5 /im). A number of works® 28 274 282 283 find qualitative agreement with the Cabrera-Mott theory, 
however linear oxide growth has been observed in other studies 263 ® 84 ® 85 *, as well as power (n < 1)E§S 
or parabolic law®^. 

Many factors can affect the experimental measurements of oxide kinetics, such as the environment 
and the type of initial copper sample. Indeed, it was found that polycrystalline film oxidation kinetics 
is almost twice as slow as single crystal oxidation kinetics 285 . The oxidation temperature is another 
factor 288 290 . O’Reilly 28 ® et al. found for a polycrystalline sample in dry synthetic air ( 02 /Ar mixture) 
that at temperatures between 250 and 500 °C the oxidation followed a cubic law, at 100 °C an inverse 
log rate and at 150 °C a linear growth rate. Gao® 9 ® et al. found a linear oxidation rate between 
200 — 250 °C with a fine-grained CU 2 O being the resultant oxide, and parabolic above that, with CuO 
the final oxidation product. 

Although a lot of work has been performed in this field, it is clear from Fig. 0 that there is 
no consistency in the results obtained. Systematic work looking at the influence of crystal structure 
(single crystal, polycrystalline, different grain size), of the experimental set-up (humidity, temperature, 
oxygen pressure) and composition of the oxide product need to be performed in order to clarify where 
the contributions to these different oxidation rates originate from. 

5.6 Native oxidation of copper in ambient conditions 

Structural details of the room temperature oxidation of copper under ambient conditions, including 
the possible influence of humidity and crystal structure, are not very clear. 

In order to study copper in ambient conditions (room-temperature, 1 atm pressure) precise measure¬ 
ment techniques are required, since surface oxide films are generally only a few nm thick. Experiments 
studying oxidation in these conditions have been performed for an enormous range of exposure times, 
between 30 minutes to 9 months, in air at ambient conditions (referred to as ambient air hereafter), 
obtaining results which are not always consistent (see Table [h]). As we will see, the outcome of the 
studies differs in the thickness of the final oxide, in the presence or not of a CuO overlayer on top of 
the CU 2 O oxide layer and in the order of formation of the oxides. The reason for this non-uniformity 
is probably to be found in the method of preparation of the film, on the type of him used (thin/bulk, 
polycrystalline/single crystal), the ambient humidity (which is rarely reported) and the surface rough¬ 
ness. 

Native oxidation has been extensively studied in copper thin films (of the order of a few hundred 
nms), because of the important role it plays in the passivation of nano-sized electronic copper parts, 
using XAS and XPS to identify the composition of the top layer and ellipsometry to measure the him 
thickness. Platzman et al. 4 proposed a three-stage oxidation mechanism involving: (a) the formation 
of a CU 2 O layer, most likely due to Cu metal ionic transport toward the oxide-oxygen interface; (b) 
the formation of a Cu(OH )2 metastable overlayer, due to the interactions of Cu ions with hydroxyl 
groups present at the surface; and (c) the transformation of the Cu(OH )2 metastable phase to a 
more stable CuO layer. Indeed, the formation of a ~ 2.0 — 5.0 CU 2 O layer first, followed by a ~ 
0.9— 1.3 CuO overlayer has been reported 3 4 291 292 . However, Lim et a/J 276 showed that the texture 
and microstructure of a thin copper him have a direct influence on the oxidation products. When 
oxidising a sample with a columnar structure and small grains the Cu 20 /CuO bylayer was obtained, 
when oxidising a uniform sample with no obvious columnar structure and grain boundaries the cuprous 
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oxide layer only was observed and oxidation was slower. 

The oxidation of bulk copper studies was studied by means of XPS, XRD and ellipsometry, and the 
results obtained are mixed. A number of studies 3 293 291 found only the formation of ~ 1.6 nm thick 
CU 2 O even after long exposure times. However, other studies of passivated Cu after exposure at room- 
temperature air 2 295 observed the formation of a few ML of CuO, which was found to start growing only 
after the CU 2 O growth process has finished. The different results obtained by these experiments can 
be due to the texture and microstructure of the copper filmEED, by sur f ace roughness or by defects 
present at the surface. Indeed, the surface roughness of metal oxide surfaces has been shown to have a 
direct influence on their wetting properties towards water, which in turn could have a direct influence 
on the formation of the native oxide 296 . In particular, very smooth copper surfaces are hydrophobic 
while rough surfaces (~ 5 nm-high roughness) are hydrophilic 297 . Furthermore the orientation of the 
crystallites at the surface influences the wetting properties and the distribution and coverage of water 
on the surface, especially in the case of copper 2 ®. Moreover, polycrystalline structures with nanosized 
grains have higher surface energy at the grain boundaries than structures that are made of micrometric 
grains or crystalline lattice, which in turns affect the wettability of the surface. 
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5.7 Nucleation sites 


A number of studies have tried to establish the role of surface defects on the nucleation of the oxide, 
in order to understand whether island nucleation is a heterogeneous process, triggered by specific 
surface features, or a homogeneous process. While some metals and semimetals can grow oxide layers 
homogeneously without the aid of impurities or surface defects (such as Ru( 0001 ) 244 , Ag( 110 ) 243 or 
Si(lll) 299 30 ®) it has been shown that other metals, like Pb 301 or Ge(lll) 3 °212Q2l j cannot oxidise without 
the presence of surface features (impurities or defects) which trigger the dissociation of O2 molecules. 

For copper it has been showd 114 161 209 255 that defect sites play a role in oxide island nucleation. 
Grain boundaries^^^, vacancy islands 209 and the edges of pits 254 i are found to be nucleating sites for is¬ 
land formation, however no preferential nucleation sites have been found at dislocations, stacking faults 
or impurities 113 254 25 ' 5 . The importance of step edges was initially inferred by Milne and Howie 2 2 ZI, 
and it was then demonstrated on Cu(lll) by means of STM and TEM 11 ' 209 . However, TEM work by 
Yang et al . 254255 showed that this is not the case on Cu( 100 ) and Cu( 110 ) films. It is indeed possi¬ 
ble, considering the different nature of oxide formation on the three low-index surfaces, that different 
defects play a more or less important role in different oxide nucleation conditions. 

More work is needed to clarify further the correlation between defect sites and oxide formation 
and remove doubts on whether small, non-structured dislocations and impurities have a role in the 
nucleation of copper oxide islands, and, from a theory point of view, why some type of defects seem 
more efficient at nucleating oxide islands than others. 

5.8 CuO formation 

CuO is expected to form after exposure of a copper surface to oxygen at high temperatures and 
pressure 39 225 . Much work has gone into understanding high temperature oxidation of copper (> 350 °) 
and it has been recently reviewed 304 305 . Most work shows the growth of CuO on top of CU2O, following 
a parabolic rate law for the thickness as a function of time. However, no studies of the atomistic details 
of the formation of CuO at high temperatures have been performed to date. 

In controlled conditions, evidence of a CuO overlayer was found by exposing a Cu sample to 
a controlled flow of O2 at high pressure and the dependence of oxidation on oxygen pressure was 
analysed. Boggio 306 investigated the pressure dependence (P = 0.03 — 7.5 Torr) on the oxidation 
of Cu(lll), and in particular the him thickness, using ellipsometry. The him growth (thought to 
be CU2O), after 90 minutes exposure at 21°, was related to the Cabrera-Mott expression of growth. 
However, a dramatic decrease in the oxidation rate with increased oxygen pressure was observed and 
related to the formation of a passivating CuO him at the oxide-oxygen interface. Pierson 226 et al. 
looked at the reactivity of a number of noble metals when subject to a how of gases. In the case of 
Cu in an O2 how, the formation of oxidised structures was found to depend on the how rate of O2. 
The formation of CU2O (at oxygen partial pressure p(C>2)=7.5 X ICE 5 Torr), CuO (p(02)=l-12 x ICE 3 
Torr) and metastable CU4O3 (p(02)=1.5 x ICE 4 Torr), was observed by XRD analysis. 

6 Conclusion and discussion 

The oxidation of copper and the physical properties of the resulting oxides are fundamental scientihc 
problems which are still not completely understood today. Since there has currently been a surge of 
interest in the use of copper oxide for catalysis, optoelectronics and gas sensing, the need for a detailed 
understanding of the surface structures of these oxides has become even more pressing. In addition, 
uncontrolled copper oxidation is still an issue in e.g. electronic applications, and understanding the 
oxide growth process is the hrst step towards mitigating it. 

In this review we have discussed the state of the knowledge regarding the structure and formation 
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of copper oxides and we have seen that the structural, optic and vibrational characteristics of bulk 
copper oxides are well understood. A good amount of computational work (mainly from DFT) has 
been performed on both CU2O and CuO, providing important information on the structures of the 
surfaces at different temperature and pressure conditions. A hexagonal structure presenting Cu surface 
vacancies was found to be the most stable on C112O and the stoichiometric (111) surface was the most 
stable for CuO. Few experimental studies are available to either confirm or disprove these suggestions 
from theory. Considering the potential technological applications of these oxides, this is an area which 
should be looked at more with experimental techniques such as STM, XPS or LEED. 

Many atomistic aspects of the formation of the oxides have been established. Indeed, the oxygen- 
induced surface reconstructions on copper surfaces with low Miller index are well known. In addition, 
STM studies have revealed certain atomic level details of the initial stages of Cu surface oxidation. 
Other experiments at low oxygen pressure show that the initial oxide growth happens through the 
formation of cuprous oxide islands which eventually coalesce. The kinetics of the oxide formation 
depends on the temperature and on the copper surface and it is dominated by oxygen surface diffusion 
and direct impingement in the first instance, and after coalescence by oxygen diffusion into the bulk. 

There is however scope for further work as many fundamental aspects of copper oxidation are still 
unclear. In fact, there is no consensus on how the transition from oxygen-reconstructed copper surfaces 
to the onset of oxidation occurs, whether this is after the surface has been reconstructed or if the islands 
start forming on the clean surface. At least for the case of Cu( 100 ) and Cu( 110 ), there is good evidence 
that the onset of oxidation will occur after the formation of an O-reconstructed overlayer which affects 
the growth kinetics of the islands themselves. The precise understanding of the formation of these 
copper islands is not only important to understand the onset of oxidation, but also to exploit the oxide 
islands as nano-templates for technological applications. In order to do that, the precise mechanism 
of nucleation needs to be understood. 

It is also unclear where the different behaviour between metals such as Ag or Al, which grow 
uniform oxide layers, and copper stems from. The O-induced surface reconstructions of different fee 
metals have been extensively studied^ 121 and the differences found in their electronic structures 307 
or surface configurations®^ could provide hints to their different behaviour upon higher exposure to 
oxygen. However, dedicated theoretical studies to this aim have not been performed so far. 

The kinetics of long-term oxide growth, both in controlled conditions and in ambient air is an issue 
that remains poorly understood. Copper oxide formation does not simply follow the Cabrera-Mott law, 
especially at the initial stages when oxide growth is not uniform across the surface. Many studies have 
been performed in order to establish the oxide formation kinetics at low temperatures, however they 
are difficult to compare to one another, because of the different conditions in use. Indeed temperature, 
humidity, oxygen partial pressure, structure of the initial copper him/coupon, the presence of defects 
and impurities seem to affect the growth of the oxide as well as the final oxide product. Moreover, the 
structure of copper oxide when grown in ambient air, which is very important to the practical uses of 
copper, is still not clearly understood, with the formation of a CuO overlayer on top of a CU2O layer 
still being debated. Further systematic experimental studies, looking at disentagling the environmental 
factors influencing oxide growth are needed. Computationally, additional studies of copper and copper 
oxide surfaces and their interaction with the environment ( e.g. water, O2, N2) would provide valuable 
information in support of the experimental work. Work in this direction has recently been reported, 
e.g. water and hydroxide adsorption on Cu and Cu oxidJ 30 ^' 312 . 

Part of the challenges that have been faced by scientists in studying these systems are due to 
the limitations of theoretical and experimental technology available to them. It is indeed clear that 
massive steps forward have been made since the first calorimetry experiments on CU2O crystals, and 
nanometre scale resolution has been obtained with STM and TEM and structural atomistic structures 
can be explored with XPS and LEED. Moreover, development of experimental techniques, such as X-ray 
lasers with extremely high spatial and temporal resolution 313 or near-ambient pressure photoelectron 
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spectroscopy 3 — might enable mechanisms of initial stages of oxidation to be explored. 

Computational techniques also present challenges from the point of view of accuracy of calculated 
structures and physical properties, as well as time length of molecular dynamics simulations. Moreover 
developments in theory, especially ab initio molecular dynamics, accelerated sampling techniques and 
more sophisticated non ab initio approaches mean that similar studies on the first steps of oxidation 
are possible. Ab initio molecular dynamics, along with a free energy sampling approach has been used 
to examine the initial stages of NaCl dissolution in liquid water 3 : it is not inconceivable that similar 
techniques could be applied to copper oxidation in ambient and even aqueous conditions. 

Finally, another important challenge, is bridging the gap between highly controlled studies per¬ 
formed in ultra-high vacuum and work aiming to understand the formation of the oxide in industrially 
relevant conditions. At the moment, these two aspects of the oxidation problem are addressed using 
different techniques at different resolutions. In order to fully understand the oxidation of copper it is 
important to be able to relate the results obtained in these different conditions to one another and 
build a unified picture of the problem. 
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